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Abstract: The aim of this article is to propose a Java concurrent program for the 
Smarandache fimction based on the equation S( p} ..... p*) = max{S(p*),...,S(p* )}. 


Some results concerning the theoretical complexity of this program are proposed. Finally, 
the experimental results of the sequential and Java programs are given in order to 
demonstrate the efficiency of the concurrent implementation. 


1. INTRODUCTION 

In this section the results used in this article are presented briefly. These concern the 
Smarandache and the main methods of its computation. The Smarandache function 
[Smarandache, 1980] is S:N*—» N defined by 


S(m) = min{k eN|kiin} (Wn EN *). (1) 
The main properties of this function are presented in the following 
(Va, be N *) (a,b) = 1=>S(a-b) = max{S(a),S(b)} (2) 
that gives us 
S(py" +...- py) = max{S(p}"),....S(pP J}. (3) 
An important inequality satisfied by the function S is 
(Va eN *) S(a)<a, the equality occurring iffa is prime. (4) 


When the number a is not prime this inequality can be improved by 


(vac N*:a not prime) S(a)< =. 


During the last few years, several implementation of The Smarandache function have been 
proposed. Ibstedt [1997, 1999] developed an algorithm based on Equation (3). The 
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implementation in U Basic provided a efficient and useful program for computing the values of S 
for large numbers. Based on it Ibstedt [1997, 1999] studied several conjectures on the 
Smarandache function. No study of the theoretical complexity has provided for this algorithm so 
far. 


The second attempt to develop a program for the Smarandache function was made by Tabirca 
[1997]. Tabirca started from Equation (1) and considered the sequence x, =k! mod n. The first 
term equal to 0 provides the value S(7). Unfortunately, the C++ implementation of this 
algorithm has been proved not to be useful because it cannot be applied for large value of n. 
Furthermore, this is not an efficient computation because the value S(7) is computed in 
O(S(n)). A study of the average complexity [Tabirca, 1997a, 1998], [Luca, 1999] gave that the 


average complexity of this algorithm is of =| : 


2. AN EFFICIENT ALGORITHM FOR THE SMARANDACHE FUNCTION 
In this section we develop an efficient version of the algorithm proposed by Ibsedt. A theoretical 
study of this algorithm is also presented. Equation (3) reduces the computation of S(7) to the 
computation of the values S(p*'),i=L...,5. The equation [Smarandache, 1980] that gives the 
value S(p* )is given by 

I p' = | i ; 

k= Yid,-—— > S(p")= 4, -P (5) 
is} P 7 i=] 


This means that if (d,,d,,,....€,) is the representation of k in the generalized base 
7] on | 
12 Pp 


ae ae i then (d,,d,,,...,d,) is the representation of S(p*) & in the generalized 
| Pp- 


i — 
base p, PP yes p’. Denote n= and 52{i]= p’ the general terms of these two bases. 


We remark that the terms of the above generalized bases satisfied: 
bil} =1, diff + J=1+4 p- diff] (6) 


bl] = p, b2fi +1] = p- bai]. (7) 
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nO 
public static long Value (final long p, final long k) { 

long 1, j, value=0; 

long b1{] = new long [1000]; long b2{] = new long [1000]; 

b1[0]=1;b2[0]=p; 

for(int 1=0:b1 [I]<=k;14+-+){b1[1+1]=1+p*b1 [1]; b2[1+ 1}=p*b2[]];} 

for(|—,j=1;j>=0;j—){d=p/b 1 fj]:p=p%b1 [j];valuet=d*b2[j]; } 

return value; 


Figure 1. Java function for S(p* ). 


Equation (5) provides an algorithm that is presented in Figure 1. At the first stage this algorithm 
finds the largest / such that 51[/] <* < bi[/ +1] and computes the generalized bases 51 and 52. 


At the second stage the algorithm determines the representation of k in the base 51 and the value 
of this representation in the base 52. 


Theorem 1. The complexity of the computation S(p*) is Ollog »P° k). 
Proof. Let us remark that the operation mumber of the function Value is 5-/, where / is the 
largest value such that b1[/] < k < b1[/ + 1]. This gives the following equivalences 


i se 
Pe egee —1sk-(p-l)<p™-le 


p-1- p-l 
op’ <k-(p-)t+i<p"™ o/silog [k-(p-I t+ </+1o 
<ol=|log, (k-(p-1)+D]. 
Therefore, the number of operations is 5- log, (&-(p —1)+1)|=Ollog, (k- p). . 


The computation of S(m) is obtained in two steps. Firstly, the prime number decomposition 
n= pi... p* is determined and all the values S(p/'),i=\...,8 are found by using a 
calling of the function Value. Secondly, the maximum computation is used to find 


max{S(p}"),...,.5(pe! )}. A complete description of this algorithm is presented in Figure 2. 
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ee Pee ee 
public static long S (final long n) { 
long d, valueMax=0, s=-1; 
if (n= =1) return 0; 
long p[] = new long [1000]; long k{] = new long [1000]; long value[] = new long [1000], 
for(d=2;d<n;d++) if (n % d = 0)f{ 
st+:p{s]=d;for(k[s]=0;n%d—=0;k[s}++,n/=d); 
value[s}=Value(p[s],k{s]); 


} 
for(§j=0;j<=s;j++) if (valueMax<value[j])valueMax=value[j]; 
return valueMax; 


Figure 2. Java function for S(7). 


Theorem 2. The complexity of the function S is Oo). 
ogn 
Proof. In order to find the prime number decomposition, all the prime numbers less than n should 


be checked. Thus, at most (7) = om) checking operations are performed [Bach & Shallit, 
: ogn ao 


1996] to find the prime divisors P>--->P, ofn. The exponents k,,...,4, of these prime numbers 
are found by k, +...+.%, divisions. An upper bound for this sum is obtained as follows 
n= p+... pi => logn=log p' -...- p* =log pt +...+logp™ = 
=k, log p, +...+k,-logp, 2k,+...+k, , 


because each logarithm is greater than 1. Thus, we have k, +...+k, <logn =O(logn). 


The computation of all the values S(pj'),j=1,...,5 gives a complexity equal to 


Slog, P,:k,;. An upper bound for this sum is provided by the following inequality 
i=l 


7§ 


log, P;*%, SA,that is true because of p,-k,<p,'". Taking the sum we find 


log 1 Pik S Yk, = O(logn), therefore the complexity of this computation is O(log 77). 
j=} 


i=} 
Finally, observe that the maximum max{S(p"*),...,S(ps! )} is determined in s<logn 


operations. 


eed n 
In conclusion, the complexity of the Smarandache function computation is oC ). > 


gn 


| 2=10000 | | 2=20000 | | n=30000 | | =40000 | | n=50000 | | n=60000 | | =70000 | | n=80000 | 


Table 1. Running times for the efficient and Tabirca’s algorithms. 


Figure 3. Graphics of the Running Times. 


Several remarks can be made after this theorem. Firstly, we have found that finding the prime 
divisors of nm represents the most expensive operation and this gives the complexity of the 


function computation. Secondly, we have obtained an algorithm with the complexity df; Z ) : 
ogn 
Therefore, this is better than the algorithm proposed by Tabirca [1988] that has the average 


complexity 2] . Table 1 shows that this algorithm also offers better running times than the 
ogn 
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algorithm proposed in [Tabirca, 1997]. These two algorithms were implemented in Java and 
executed on PENTIUM II machine. The times [milliseconds] of the computation for all the values 
S(i) , i=1,...,.n were found, where n=10000,....,80000. Row 4 1 gives the times for this efficient 
algorithm and row A 2 gives the times for the algorithm proposed in [Tabirca, 1999]. Another 
important remark drawn from Table 1 is that the difference between the times of each column 
does not increase faster [see Figure 3]. This is happen because the complexity of the algorithm 
proposed by Tabirca [1997] is of =) : 


3. JAVA CONCURRENT ALGORITHM FOR THE SMARANDACHE FUNCTION 

In this section we present a Java concurrent program for the computation described in Section 2. 
Firstly, remark that many operations of this algorithm can be performed in parallel. Consider that 
we know all the prime numbers less than 7. Usually, this can be done by using special libraries. 


Let D,,....P, be these numbers. Therefore, we can concurrently execute the computation of the 


exponent of p, and the computation of the value S(p* ). 


A Java program may contain sections of code that are executed simultaneously. An independent 
section of code is known as a thread or lightweight process [Smith, 1999]. The implementation 
presented here is based on equation (3): S(pj'-...-p** )=max{S(pi*),...,S(p® )}. Each 
S(p') is calculated concurrently in a thread. On single processor systems, the use of threads 
simulates the concurrent execution of some piece of sequential code. The worst case execution 
time can be taken as the longest execution time for a single thread. On a multi processor system, 
given enough processors, each thread should ideally be allocated to a processor. If there are not 
enough processors available, threads will be allocated to processors in groups. Unlike pure 
concurrent processes, threads are used to simulate concurrency within a single program. Most 
current everyday programs use threads to handle different tasks. When we click a save icon on a 
word processing document typically a thread is created to handle the actual saving action. This 
allows the user to continue working on the document while another process (thread in this case) is 
writing the file to disk. 


For the concurrent algorithm consider the Java function for S(7) in Figure 2. Typical areas that 


can be executed concurrently can be found in many loops, where successive iterations of the loop 
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do not depend on results of previous iterations. In Figure 4, we adapt the for loop (Figure 2) to 
execute the Value function (Figure 1), responsible for calculating S (p*), concurrently by 
creating and executing a ValueThread object. When all the required threads have begun 
execution, the value of max will not be known until they have completed. To detect this, a simple 
counter mechanism is employed. As threads are created the counter is incremented and as threads 
complete their tasks the counter is decremented. All threads are completed when this counter 


reaches 0. 


public long S(iong n) 
af 
if (n==1) return (long); 


Prime decom = new Prime(n); 
noPrimes=decom.noPrime(); 
if (noPrimes = 0) 

value = null; 
value = new long[noPrimes]; 


for (int k=0;k<noPrimes;k++) 
{ 
started; 
new ValueThread(decom.getPrime(k), decom.getPow(k), this, k); 
} 
while (started > 0) 
{ 
c 
Thread. yield); 
}catch (Exception e) 
{ 
} 
} 
return max; 


Figure 4. Modified Java function for S(7), used to concurrently execute the Value function 


As each. thread completes its task it executes a callback method, addValue (Figure 5). This 
method is declared as synchronized to prevent multiple threads calling the addValue method at 
the same time. Should this be allowed to occur, an incorrect value of the number of threads 
executing would be created. Execution of this method causes the value array declared in method 
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S (Figure 3) to be filled. This value array will only be completely filled after the last thread makes 
a call to the add Value method. At this point, the value of max can be determined. 


public synchronized void addValue(int k, long val) 
{ 


for (int i=1; i<=k; i++) 
if (value[i] > max) 
max = value[i]; 


3 . 3 
a ES ee et 


Figure 5. The addValue method called by a Thread when its task is completed. 


This algorithm illustrates how concurrency can be employed to improve execution time. It is also 
possible to parallelise the algorithm at a higher level, by executing the function responsible for 
calculating each S(7) in an independent thread also. Tests of this mechanism however show that 


it is more efficient to only parallelise the execution of 5. (p*). 


The concurrent Java program has been run on a SGI Origin 2000 parallel machine with 16 
processors. The execution was done with 1, 2, 4 processors only and the execution times are 
shown in Table 1. The first line of Table 1 shows the running times for Algorithm Al on this 
machine. The next three lines present the running times for the concurrent Java program when 
p=1, p=3 and p=4 processors are used. 


| n=20000 | | n=30000 | 0000 | n=40000 | n=50000 | 2=60000 | | n=70000 | | 2=80000 0000 
Al 9832 19703 31237 49774 68414 96242 115679 
Em Se a a aS 
CA (p=1) 9721 19474 30195 49412 68072 95727 115161 
CA (p=2) 5786 11238 22872 31928 42825 60326 75659 
ee) GU eas | sive ees | fees 
CA (p=4) 3863 7881 14017 19150 30731 42508 53817 


Table 2. Running Times for the Concurrent Program. 
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4. CONCLUSSIONS 
Several remarks can be drawn after this study. Firstly, Equation (3) represents the source of any 
efficient implementation of the Smarandache function. In Section 2 we have proposed a 


sequential algorithm with the complexity f=.) . We have also proved both theoretically and 
ogn 


practically that this algorithm is better that the algorithm developed in [Tabirca, 1997]. 


Secondly, we have developed a Java concurrent program in order to decrease the computation 
time. Based on the thread technique we have performed concurrently the computation of the 
values S(p*). This concurrent implementation has proved to be better than the sequential one. 


Even running with one single processor the times of the concurrent Java program were found 
better than the times of the sequential program. 
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Appendix A 


The full code for the concurrent implementation presented in Section 3. 


// Smarandache.java 


import java.io.*; | 
import java.util. *; 


public class Smarandache 


{ 
public Smarandache() 
{ 


long n=0, L, j; 
long val; 
BufferedReader br = new BufferedReader(new InputStreamReader(System.in)); 
try . 
{ 
System.out.print ("n ="); 


n = Integer.parseInt(br.readLine()); 
}catch (IOException e) 
{ 


System.out.printin ("IOException : "+e.getMessage()); 
System.exit(1); 
} 


Smar sm = new Smar(); 


Date begin = new Date(); 
for (i=1; i<=n; i++) 


{ 


} 
Date end = new Date(); 
System.out.printin ("Time good is "+ (end.getTime() - begin. getTime())); 


val = sm.S(i); 
} 


public static void main (String args[]) 
{ 


new Smarandache(); 
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/{ Smar java 


public class Smar 


{ 


private long value[]; 

private long max = Long. MIN_VALUE; 
private int noPrimes=0; 

private int started = 0; 


public Smar{) 


public long S(long n) 


if (n==1) 
return (long) 0; 


Prime decom = new Prime(n); 
noPrimes=decom.noPrime(); 


if (noPrimes = 0) 

value = null; 
value = new long[noPrimes]; 
for (int k=0;k<noPrimes;k++) 
{ 

started; 

new ValueThread(decom. getPrime(k), decom. getPow/(k), this, k); 
} 
while (started > 0) 
{ 

c 

Thread. yieldQ; 

}catch (Exception e) 

{ 

} 
} 
return max; 
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public synchronized void addValue(int k, long val) 


{ 
value[k] = val; 
started--; 
if (started == 0) 
{ 
max = value[0]; 
for (int i=1; i<=k; i++) 
if (value[i] > max) 
max = value[i]; 
} 
} 
} 
a eee rs es ee ee ee 
//Prime java 
public class Prime 
{ 
private int s; 


private long p[]=new long [1000]; 
private int ord[]=new int [1000]; 


public Prime() 
{ 
s—0; 
} 
public Prime(long n) 
{ 
long d; 
for(d=2,s=0;d<=n;d4+) 
if(n%d = 0) 
{ 
p{s}=d; 
for(ord[s]=0;;ord[s}]++,n=n/d) {if(n%d!=0)break;}; 
st; 
} 
} 
public int noPrime() 
{ 
return S; 
} 
public long getPrimetint i) 
{ 
return p[i]; 
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public int getPow(int i) 
{ 


return ord{i]; 
} 
} 
// ValueThread java 
public class ValueThread 
{ 
private long p=0, a=0; 
private Smar owner; 
private int index = 0; 
public ValueThread (long p, long a, Smar owner, int index) 
{ 
this.p = p; 
this.a = a; 
this.owner = owner; 
this.index = index; 
run(); 
} 
public long pseuPow(long p, long a) 
{ 
if (a = 1) 
return (long) 1; 
return 1+p*pseuPow(p,a-1); 
public long Pow(long p, long a) 
{ 
if (a == 1) 
return (long) p; . 
retum p*Pow(p,a-1); 
} 2 
public void run() 
{ 
long rest=a, val=0; 
int k, i; 
for(k=1 ;pseuPow(p,k)<=a;k++);k-; 
for(i=k;i>0;i—-) 
{ 
val += Pow(p,i)* (long)(rest / pseuPow/(p,i)); 
rest %= pseuPow(p,i); 
} 
owner.add Value(index, val); 
} 
} 
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